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1 Introduction 

In the past few years, first-passage properties in nonequilibrium situations have attracted 
much interest, with numerous works devoted to the so-called persistence phenomenon (for 
a recent review, see |]I|, and references therein). The simplest setting might be within 
the field of phase-ordering kinetics where, as domain growth proceeds, the following 
question naturally arises 0, |^: What is the fraction R(t) of space which has always been 
in the same phase (say the + one) up to time t ? Equivalently for zero-temperature 
quenches, R(t) is the probability that a given spin has never fiipped up to time t, and is 
observed to decay in time as t^^ , 9 being the persistence exponent. The latter can also be 
considered as a first-passage exponent, since — di?(t)/dt represents the probability that 
an interface first sweeps over a given point at time t. Experimental measurements of 9 
have been conducted in breath figures |^, and in a system of nematic liquid crystal akin 
to the 2d Ising model dynamics [Q. At par with the general view that thermal noise 
should be asymptotically irrelevant in coarsening systems, several numerical methods 
P) i) S Hi indicate that 9 seems to keep the same value for all quenches in the broken- 
symmetry region. Unfortunately, exact calculations of this new critical exponent are 
scarce |^, |10[ since, after integrating out at a given location in space the spatial degrees 
of freedom, one has to deal with an effective stochastic process which is generically 
non-Markovian. Perturbative techniques to calculate 9 have therefore been developed 



TT| , |T^, [T^, |lj, |T5], especially when a spin can be associated with (the sign of) a Gaussian 



process, a situation common to various closure schemes ||16| , \u[ ^8| of phase-ordering. 
Other lines of investigation have been concerned with extending the notion of (local) 
persistence probability. This includes notably the definition of a persistence exponent 9g 
for the global order parameter (at Tc [|l^, and below @, |^), the notion of block-persistence 
0, 01 (which encompasses both 9 and 9g), and two others generalizations [^ ^, relaxing 
in a different fashion the no spin-flip constraint inherent to the deflnition of R{t). Here 
we shall deal with the extension introduced in [^ (see also 1^), which is based on the 
consideration of the local mean magnetization: 

Mit) = - f&t'ait'). (1.1) 

t Jo 

This quantity is simply related to the proportion of time spent in one the two possible 
phases by a given spin a = Oy.. The distribution of M(t) contains also the persistence 
probability, for Prob{M(t) = +1} = i?(t). Yet, as soon as |M(t)| < 1, changes of phase 
are allowed, and one observes that M(t) converges, when t —>■ oo, to a broad limiting 
distribution Jm, reminiscent of the classical arcsine law for Brownian motion, whose 
singular behaviour Jm^x) ~ (1 — x^Y^^ at the edges x = ±1 of its support involves the 



persistence exponent ||2^. Another conceptually related route giving access to 9 is to 
consider, for a fixed value of the level x < 1, the persistent large deviations R{t,x) = 
Prob {M{t') > X, Vt' < t} of the mean magnetization, which can be also interpreted as a 
the (+)-persistence of the process sign[M(t) —x]. For large times, this quantity decays as 
a power-law t~^^^^ pO|, and the original persistence exponent appears as the limiting case 
6 = linia;^!- 6{x) of a continuously-varying family of generalized persistence exponents. 
Let us mention that another infinite hierarchy of exponents appears when one allows 



a spin to be "reborn" with a certain probability p each time it flips pTf, or when one 



considers the survival properties of (spatial) domains in coarsening systems [Q. In 
fact, R{t, x) displays a non-trivial structure even for the simple (uncorrelated) binomial 
random walk 124 



Although the existence of level-dependent exponents is by no means new (see, e.g.. 



25| for Gaussian processes, and the more recent work [|T5[), it is physically very appealing 
that, as pointed out in pf, both the limit law fui^x) and the spectra 0{x) seems to be 
universal for the 2d Ising model in the whole region T < T^., up to a rescaling of x 
by the corresponding equilibrium value of the magnetization, thus reflecting the fact 
that coarsening systems, though perpetually out-of-equilibrium, are nevertheless in local 
equilibrium. 

We have also argued recently [^ that the global shape of the 6{x) curve is a sensitive 
probe of the type of noise (stochastic or deterministic) vehicled by the domain walls' mo- 
tion, thereby suggesting an explanation to the discrepancy currently reported 0] between 
the persistence exponents of the 2d Ising model and of its continuous counterpart, the 
time-dependent Ginzburg-Landau (TDGL) equation. 

The purpose of the present work is to present a family of models for which gener- 
alized persistence properties can be obtained analytically. We shall do this within the 



framework of the independent interval approximation (IIA) used in [|T2i [T3| to calculate 



the persistence exponent for a Gaussian random field evolving under the diffusion equa- 



tion. The IIA has been used later to determine the limit law for this system ||20|, and 
our work is a natural continuation of this study. Although the IIA is an uncontrolled 
approximation, limited to "smooth" processes (that is, if a{t) = signX(t), the process 
{X{t)}^ must be everywhere (with respect to t) differentiable), this method has given 
very accurate results, both for 6 and fM{x). Excellent agreement will be also found here 
when comparing our results with the spectrum 6{x) obtained numerically for the diffu- 
sion equation. In a certain limit, there will also appear an unexpected similarity with 
an exact relationship obtained for the d{x) exponents of the Levy-based model studied 



in [^, although the latter process resembles the Id Ising model and is non-smooth in 



nature. To conclude our motivation, we mention that some of our analytical results 
rely on a technique to deal with multiplicative random recursions (defining a so-called 
Kesten variable [EH]), on which constraints are enforced. Since such variables show up in 



a variety of one-dimensional random systems [^ ^, it is also hoped that the method 
introduced here might prove helpful to tackle new questions in this field. 

The paper is organized as follows. In Section 2, we review the IIA. Section 3 and 
Section 4 are respectively concerned with the determination of generalized persistence 
properties when the number of spin-fiips is fixed (Section 3), and when the time is fixed 
(Section 4). We conclude in Section 5 with a discussion. 

2 The independent interval approximation 

To fix notations, we first give a brief account of this method used in [|l^, [l^ to calculate 
the persistence exponent for the diffusion equation. Owing to dynamic scaling, coars- 
ening systems with algebraic growth laws can be rendered stationary by going over to 
"logarithmic" time £ = In t — Into- The principle of the IIA is to assume that the lengths 
of the (logarithmic) times /j = Intj — lntj__i between successive spin-fiips ti,t2, . . . form 
a renewal process, the common probability distribution function /(/) of these intervals 
being determined by relating it to the scaling form {a(to)a(t)) = a(toA) = A{£) of the 
spin-spin autocorrelation function (which is known exactly for the diffusion equation, 
and more generally for any Gaussian process). The relationship reads, in Laplace space 
(/>)=Cd£e-V(^)), 

/(.) ^ ^-(^)i^-^Ms)]/2_ 
1 + (0 [1 - sA{s)]/2 
In ( |2.1|) , the average length (/) between two zero-crossings is determined by the small-£ 
behavior A{i) ^ 1 ~ 2i/ (/) of the correlator, li A{i) — 1 ~ i° with an exponent a < 1, as 
happens in the Id Ising model (for which a = 1/2) or for some interface growth models 
|l3T| , ^ (where a is related to the roughness exponent), then A'{0) is infinite and the IIA 
cannot be used as such (see, however, ||33[) to enumerate the spin-fiips: any change of 
sign is likely to be followed by a dense sequence of zero-crossings (of fractal dimension 
Do = 1 — a IQ). Otherwise, the expected large-time behavior R(t) ~ t~^ oc e~^^ of the 
persistence probability can be inferred from the existence of a pole at s = —6 for f{s) 
||35|| . This procedure can be implemented in the Ohta-Jasnow-Kawasaki closure scheme 
of phase ordering |jl6| (abbreviated from now on as the "diffusion equation"), where the 
dynamics of a spin is recast in terms of an auxiliary diffusing field, 

^0(r,t) = DVV(r,t), (2.2) 



evolving from zero-mean random (Gaussian, or simply short-ranged) initial conditions. 
The zeros of this diffusing field are interpreted as representing the positions of do- 
main walls: o"r(t) = sign0(r, t). At a given location in space, the normalized process 
X(t) = (/)(r, t)/[(0^(r, t))]^/^ is still Gaussian, therefore solely determined by its two-time 
correlation function: 

-|d/4 



c(to,t) = (X(to)X(t)) 



4tnt 



(2.3) 



which depends only on the ratio to/t of the two times. The spin-spin autocorrelation 
function is ((T(to)o"(i)) = (S/vr) arcsinc(to,^), (this equality holding true for any Gaussian 
process), and reads, when expressed in terms of the stationary timescale £ = In t — Into, 

A{£) = -arcsin|[cosh(£/2)]-'^/2\ ^2.4) 

Plugging this back into (|2.1|) gives values of the persistence exponent 6 extremely close 



compared with those obtained by a direct simulation of (|2.2|) : one finds respectively 



in rf = 1,2,3 space dimensions Sua = 0.1203, 0.1862, 0.2358, while Onum = 0.1207(5), 



0.1875(10), 0.2380(15) |T2|, |T3 



The IIA can also be used to calculate the limit law of the mean magnetization [pO 
For our later purposes, it is convenient to introduce the ratios Zi = ti_i/ti = e~'' of two 
successive flip-times. We denote their common distribution by fz- The Laplace trans- 
form f{s) of the probability distribution of the /j's is nothing but the Mellin transform 
jQdzz^^^fz{z) of that of the Zi^s, which implies that 

fziz)^aoz'-\ z^O. (2.5) 

In ( |2.5| ), the proportionality constant Oq can be computed from the residue at s = —6 of 
the Mellin-Laplace transform ( |2.1| )). According to the sign of o"(t), one can now rewrite 
tM{t) = ±[ti- (t2 - ti) + . . .] as 

M{t) = ± (l - 2e-^Wn) . (2.6) 

Here, ^ = In t — In t„ is the length of (logarithmic) time since the last zero-crossing, and 
the "weights" WnS (which belong to the interval (0, 1)) are identified as Kesten variables, 
for they obey the following random multiplicative relation: 

Wn = l-ZnWn-l. (2.7) 

In the large-time limit, the random variable Wn converges in distribution to the solution 
oi W = 1 — ZW, which means that the probability distribution g oi W has to satisfy the 
following integral equation: 

g{w) = C dw'w'-'g{w')fz (^^) . (2.8) 



The solution of ( |2.8| ) represents the major step for the determination of the hmit law of 
the mean magnetization, since for a renewal process the law of the backward recurrence 
time ^ is known: (e~^^ ) = (1 — /(s))/ (/) s. One of the very few cases where the Dyson- 
Schmidt integral equation (|2.8|) can be inverted explicitly is when the distribution fz is 
a pure powerlaw in z, that is when fz{z) = 6z^~^, \fz G (0, 1) |3^. One then finds 



g{w) = B-\e,e + i){i-wy-'w\ (2.9) 

{B{a,b) = r(a)r(6)/r(a + b) being the usual Beta function), and eventually that the 
limit density of the mean magnetization is a Beta-law on (—1, 1): 

fM{x)=B-\l/2,9){l-x'Y-\ (2.10) 

This particular choice for fz, proportional to the tail-behavior of the true distribution, 
nevertheless contains much of the physics of the persistence phenomenon since, in all 
coarsening systems, the limit density of the mean magnetization is found to be numeri- 
cally indistinguishable on its whole support from the law ( p.lO|) pO| , |2^ , H, ^ . This can 



also be demonstrated within the llA for the diffusion equation (using the full law fz{z) 



determined by (|2.1|) - (2^)), where one finds (by extrapolating the large-order behavior of 
the moments of M) that the tail-behavior JmIx) ~ (1 — x^)^^^ contains an overwhelming 



proportion of the full measure [ 20 



Let us now define the quantities of interest in the present work. Depending whether 
the number of spin-flips n or the time t is fixed, we want to calculate 

R{n,x) = ProbjMfc >x, V A; = l,2,...,n}, (2.11) 

R{t,x) = Prob {M(t') > X, V t' < t} . (2.12) 

In ( p.ll| ), Mk = M{tk) is the value of the mean magnetization at exactly the k-th flip- 
time. These two aspects of generalized persistence properties are different, but we shall 
show that they are related. 

3 Persistent large deviations of the mean magneti- 
zation when the number of spin-flips is flxed 

3.1 General formalism 

Let us start with a study of R{n, x). We assume from now on that initially cr(0) = +1, 
and that the the first spin- flip ratio zi possesses the same law as the others: zi = to/ti, 
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with to = 1- Since the last condition {M„ > x} in ( p. 11 ) can only be violated — had it 



been satisfied earlier — when a spin is journeying in the minus phase, that is, after an 
odd number of spin-fiips, one has R{2n + l,x) = R{2n,x). Expressed in terms of the 
weights W2k introduced earlier, the n conditions to be maintained at the 2nth spin-fiip 
time are: 

R{2n, x) = Prob |w^2 < ^^, W^ < ^^, . . . ,W2n < ^y^} • (3.1) 

Instead of considering all the instances taken by the W2kS and then applying the con- 
straints (^.11), it is sufficient to consider the ensemble of constrained variables ^2^ G 
[0, i^] transformed by ^2^ = 1 — Z2A:(1 — ^2/0-1^2^-2) (henceforth to simplify notations 
we do not display the dependence of the variables W2k upon x). This does not define 
a map on [0, ^^], since some variables may fall outside the interval. Let us introduce 
P2kix), which represents the fraction of the variables which are left on [0, ^-^] during 
iteration 2k. Therefore one has, 

n 

Ri2n,x) = l[p2kix). (3.2) 

k=l 

The key of the argument is now to write the dynamics for the distributions of the 
constrained variables W2k on the interval [0,^-^]. Their renormalized densities h2k 

l-x 

(/o ^ dwh2k{w) = 1) have now to obey 

Jl-w W' \ W' / Jl-w' W" \ W' J 

The subsequent step is to recognize that in the large-ra limit the renormalized densities 
/i2n's have to reach a stationary distribution h{w), vanishing identically for w > ^^. 
Once this stationary regime is reached, one has simply 

i?(2n,x)~[p(x)]", n>l, (3.4) 

with p{x) = lim„^ooP2n(a^) determined self-consistently through the solution of 

p(.)ft(,„) = /' ^4f, (l^) f-^^^hWV. (i^) . (3.5) 

Jl-w W \ W' J Jl-w' W \ W' J 

subject to the constraint that h{w) be a normalized probability distribution on the in- 
terval [0, ^-^]- We have to stress that this exponential decay of R{n,x) is in agreement 
with the expected algebraic behavior of R{t, x) since, for a smooth process, the typical 
number of spin-fiips up to time t scales as In t/ (/) . This would not be true anymore 
for the l(i-Ising model (or for the Levy-based model studied in ||2^), where the typical 
number of spin-fiips rather behaves as t^'^ |3^ (or t^ for the Levy model), and R{n,x) 
instead goes to zero algebraically. 



3.2 Exact solution when fz{z) = 9z^ ^ 

Since the determination of the unconstrained probabihty density g of the W^s is aheady 
a difficult problem, there is no hope to solve ( |3.5| ) for an arbitrary distribution fz- Never- 
theless, when fz{z) = 9z^~^, the structure of ( p.5|) still allows for an exact calculation of 
p{x), which we now present. [Note that for the probability distribution function /(/) of 
intervals between spin-flips on the logarithmic timescale, this choice for fz{z) corresponds 
to a Poisson process: /(/) = 9e~^^.] In this case, rewriting ( |3.5|) as 

two successive differentiations with respect to w show that h{w) obeys the following 
(purely local) differential equation 

d^h dh r 1 

w{l-w)^ + [l-9-{3- 29)w] 3- - (^ - 1)' - 9yp{x)] h = 0, (3.7) 

aw^ aw L J 

which is recognized as Gauss' hypergeometric differential equation. The integral equation 
(TE) demands that h{0) = 0, and this selects the solution 



9 9 
h{w) = NoTm.w\Fi \ 1 ^=, 1 + ; l + 9;w\ , (3.8) 

where 2-^1(0, b; c; w) is the hypergeometric series: 

r, ,, , _ r(c) ^ r(a + n)r(b + n) w" 

^^''''•''■'^■"'' = fMrWS r(c+n) ^- P-»' 

[One remarks that for x = —1, p(— 1) = 1, and (|3.8| ) does reduce to ( |2.9| ), since 
2-Fi(l — 9,1 + 9;1 + 9;w) = (1 — wY^^.] To determine p{x), we need a second boundary 
condition, which is obtained by differentiating the left-hand-side of ( |3.6| ) with respect 
to w, whereupon we set w = ^-^. Then the right-hand-side of the obtained equation is 
identically zero, since the range of the remaining integration vanishes. Using the formula 



d 

dw 



a+fe+c-l„,,c-2 



{l-w)''+''-''w''-'2Fi{a,b;c;w) = -{I - c){l - w)''+'+''-'w'-' x 



we obtain, for ^ 7^ 1, 



X 2Fi{a-l,b-l]c-l]w), (3.10) 



I —9 9 „ 1 — X , , , 

2F,\^=,^=;9; -^— \ =0. (3.11) 



.\/p(aO' \fp{x) 



2 




Figure 1: The decay rate p{x) of R{2n,x) (as defined by (|3.4D) for three representative 
values of 6': 6 = 0.121 (the value of the persistence exponent for the Id diffusion equation), 
e = 0.2 (that of the 2d TDGL equation [|, ||), and ^ = 0.5, from bottom to top. Full 



line: exact formula (p.ll|), circles: numerical simulations. 



[For 6 = 1, Eq. ( |3.10| ) is not valid, but a similar identity applies, yielding 2-^1(1 — 



p~'/^{x) 



1+P 



-1/2, 



X) 



2; ^2^) = 0.] Equation ( p.ll| ) determines implicitly p{x) as a func- 



tion of X, for a given value of 6. This formula can be inverted when 6 = 1/2, to give 

(3.12) 



p(x) 




{6 = 1/2). 



In this case, the constrained density h{w) has also an explicit expression: 



2 arcsin J^-^ 
h{w) = Norm. , — sm 



vr arcsm ww 



TT 



VT 



w 



2 arcsin ji^^ 



{9 = 1/2), 



(3.13) 



where the normalization factor (which depends of course on x) can also be calculated. For 
other values of 6, ( |3.11| ) is generically multivalued, but there is a unique determination 
of p{x) in the physical range (0, 1), and the ensuing solution can be obtained numerically 
with an arbitrary precision. 

Figure |I| shows our analytical prediction ( p.llD for three different values of ^, compared 
with the results of a direct numerical simulation of the process defined by (p.7|)-(pTTD, the 
ZnS being distributed according to fz{z) = 9z^~^. The agreement found is excellent. 
Numerical measurements of the expected exponential decay of R{n, x) for the diffusion 
equation are hard to perform, for the number of spin- flips scales as Int/ (/), which is only 
of the order or 2 — 3 even for times t = 10^. 



Even if (|3.11|) is rather unwieldy for a generic value of 9, one can extract from it the 
behavior of p{x) at the two edges x = ±1 of the spectrum. This turns out to be also 
possible for a generic law fz-, and we shall now sketch these calculations. 

3.3 Limiting behaviour of p{x) when x ^ 1 



When X ^ 1, p{x) -^ 0, and ( 3.11 )-(3.9) show that the proper scaling variable is S* = 
(1 — x)/2p{x). In this scaling limit, we have to keep all terms in the hypergeometric 
series ( |3.8| ), which simplifies to 

0-Ep(^^-mS<-VV._.(2v^). (3.14) 



where Je-i{v) is a Bessel function. The solution of (|3.14|) is therefore 



, , ( 29 V l-x , , 

P(^) ^[- ^, (3.15) 

where j6»-i,i is the first strictly positive zero of the Bessel function Jq^i{v). [One can 
check with (|3.12|) that this result is correct for 9 = 1/2, since in this case one has simply 



J-i/2{v) = cosf.] As expected, the decay rate — Inp(x) of the generalized persistence 
probability R{n,x) diverges when x — >■ 1, since no spin-flips are allowed in this limit. 

It is also possible to estimate the behavior of p{x) when x —>■ 1 for the "true" law fz 
given by the inverse transform of (|2.1|) . To do this, we need to determine the behavior 



of fz{z) near z = 1. This can be done by developping (|2.1|) for s — >■ oo in powers of 1/s. 
The resulting expression involves the coefficients of the Taylor series expansion of the 
field-correlator near i = 0: C{i) = {X{0)X{i)) = 1 + 02^72! + £4^74! + . . .. To lowest 
order in 2; — ;> 1, the result for the diffusion equation in d dimensions is 

fz{z) ^ ai(l - z), with a, = ^^ = ^. (3.16) 

A careful study of the integral equation ( p.5|) , using the fact that the singular part of 
g{w) in the limits w — ;> 0, 1 is still given by { \i.9\) for any law fz{z) ~ aoZ^~^, shows 
eventually that p{x) ~ (1 — x)^ when x ^ 1. 



3.4 Limiting behaviour of p{x) when x ^ —1 

In the opposite limit x ^ —1, p{x) -^ 1, and one can derive, by retaining the first 



two terms in the series ( |3.11D (the natural expansion parameters being 1 — J p{x) and 

10 



y = (1 + x)/2), the following expression: 



1 - dp{x) ^ y" 



ev^ie) e^T\e) 



y~r o\T- o\ ~ y (1 + ^^ ^O^ ^^) 



T{2e) ^{i-e)v{2e) 



-1 



(3.17) 



]-B-\e,e + i)(^^) (3.18) 



[The higher-order correction term in ( p.l7| ), even though it looks superficially divergent 



for 9 = 1, can be analytically continued to give the correct result 



l+x 

p{x)^ ^ 7 T. — ^, ^^1. ^=1-] (3-19) 



2 2^2 

To lowest non-trivial order in (l + x)/2, using (|2.9|) , ( p.l8| ) is therefore equivalent to 



/ 1 I \ ^ 1 — a: 

p(a;) ^1-5-^(^,^ + 1) (^^j ^ f ' dw g{w), x -^ -1. (3.20) 

This result has a simple and intuitive interpretation: in this weakly-constrained limit, the 
W2kS obeying ( |3.1D are simply asymptotically distributed according to the unconstrained 
probability distribution g{w). 

Again, such a behavior is not restricted to the particular choice fz{z) = 9z^^^: for 
an arbitrary distribution behaving as in ( |2.5| ), one can show, still using the fact that 
g{w) ~ w^{l — wY^^ when w — > 0, 1, that the Ansatz 

l-p(x)~(i±^)', x^-l, (3.21) 

does solve the integral equation (|3.5|) . 



4 Persistent large deviations of the mean magneti- 
zation when the time is fixed 

Even for the simplest choice fz{z) = 9z^~^ (corresponding to a Poisson process with 
mean 1/9 on the logarithmic timescale), we have not succeeded in calculating exactly the 
complete spectra 9{x) of generalized persistence exponents: when expressed in terms of 
the spin-fiip ratios, the constraints that one has to enforce are non-linear and non-local. 
For an arbitrary law fz-, we can nevertheless offer a systematic perturbative development 
in the neighbourhood of a; = 1, and a non-perturbative resummation in the vicinity of 
X = — 1, the latter being presumably exact (though we cannot strictly prove it). 
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4.1 Limiting behavior of 9{x) when x ^ 1 

The first idea is to follow by continuity 6{x) from 6 = 6{1) in an expansion in the number 
of spin-flips which have occurred up to time t. The first term of the expansion for R{t, x) 
is trivial, and just corresponds to the persistence probability. That is, if no spin-flips 
have occurred up to time t, then the first flip-time ti = 1/zi > t, and the probability 
of this event is Rq = Jq dzifz{zi). If just one spin-flip takes place up to time t, then 
the two first flip-times satisfy ti < t < t25 or equivalently the two first flip-ratios obey 
t-^ < zi < 1, < Z2 < (tzi)-^. The magnetization M(t) = f-^ (ti - (t - ti)] = -l + 2ti/t 
has to satisfy M{t) > x to contribute to R{t,x). This tightens the upper bound on the 
first flip-time ratio zi, and gives a contribution 

/2i = / dzJziz^) / dz2fz{z2) (4.1) 

Ji/t Jo 

to R{t,x) (as in Section 3, we still denote y = [1 + x)/2; note also that the IIA has 
been made by writing the joint probability distribution of Zi,Z2 in factorized form). 
Consider now the situation where a spin has experienced two spin-flips up to time t. 
This corresponds to the following bounds on Zi,Z2,Zs: 

r^ < Zl < 1, (tzi)-' <Z2<1, 0<Z3< {tziZ2)-\ (4.2) 

There are now two contributions to R{t, x) to distinguish, according to the range of values 
taken by the first flip-time ratio zi. If the latter obeys the same bounds as in ( [4.1D , then, 
irrespectively of the subsequent flip-time sequence ^2, ^3, it is not hard to check that such 
an instance will always contribute to R{t,x), and the bounds on Z2, z^ are just those 
defined in (|3): 

rl/{yt) rl r{tziZ2)-^ 

R2.1 = / dz^fziz^) / dz2fziz2) / dz^fzizs) (4.3) 

Jl/t J(tz-i_)-'^ Jo 



"l/iyt) rl r{tz^Z2)-^ 

dzifz{zi) / ( 

'l/t J{tzi)-^ 

The remaining contribution to R{t, x) after two spin-flips reads 

"1 /•! /•(t^i^2)-l 

dzjzizi) 
'i/{yt) Jy 



R2.2 = t dzifz{zi) f'dz2fz{z2) I *'''' dz^fziz:,). (4.4) 

Jl/ivt) Jy Jo 



We now evaluate the sum of these four contributions in the long-time limit, and for 
y = (1 + x)/2 close to 1. The derivation will show that we do not need to take into 
account higher-order contributions in the number of spin-flips if we content ourselves 
with an expansion of 6{x) — 9 to lowest non-trivial order in (1 — x)/2 = 1 — y ^ 0. In 



most of the integrals, the large-t behavior is determined by the small- 2; form (|2.5|) of the 



distribution function of the ratio of flip-times, while other spin-flip ratios can be summed 
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Ri ^ ao r^'^z^z,''^ I'dz2fz{z2) ^ ^{y-' - 1) ^ ^(1 - y). (4.5) 



over without restriction due to the fact that y ^ 1. For instance, to lowest-order, ( |4.1| ) 
simphfies to 

/t 

This kind of simphfication is also at work for Rq and R2,i, which are also found to be 
simply proportional to t'^ with y-dependent amplitudes, but not for -R2,2, where the 
intermediate integral in equation (^4.4|) shows that what matters is also the behavior 
( 3.16| ) of fz{z) when z —>■ 1. One eventually finds that 

R.. - ^^i^#^. (4.6) 

The appearance of a term involving a logarithm of the time is mandatory to determine 
9{x), for if we develop the expected algebraic behavior R(t,x) ~ t~^(^) for both t ^ 1 
and X — > 1, we would obtain 

f-e{x) _ exp[-^(l)lnt+(l-a;)^'(l)lnt + ...] 

^ r^[l-{l-x)e'{l)lnt + ...]. (4.7) 

Comparing the two expansions (|4.6|) and ([4.7|), one finds that for x — *> 1 



e{x) 



a^ai ^1 — x^ ^ 



(4. 



20 V 2 
where Oq and ai are the coefficients describing the behavior of fz{z) respectively near 



z = and z = 1 (equations (p.5| ) and (|3.16|) ) 



4.2 Limiting behavior of 9{x) when x ^ —1 

The preceding subsection has shown that there is an interplay between the distributions 
of the flip-time ratios, the level x, and the time t, which appears hard to handle in a 
systematic fashion. More precisely, the distribution of the number N^ of spin-flips which 
have occurred up to time t and which contribute to the generalized persistence probability 
R{t, x), depends on the level x. A scaling argument such as (N^) ~ ( A^^ \ = Int/ (/) 
— although qualititavely correct to explain the exponential dependence of R{n, x) and 
the converse algebraic one of R{t, x) — does not suffice for a quantitative calculation. 
However, we expect that for x ^ —1 we can use the unconstrained distribution for the 
number Nt = NJ; of spin-flips which have taken place up to time t, providing that we 
take into account the correct form of the survival probability R{2n, x) every other time 
step. Namely, one should have the decoupling equation: 

Prob {M(t') > X, \ft' <t} ^ Prob {M(tfc) > x, Wk < n} x Prob {Nt = n} . (4.9) 
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Taking the equation (^4.9|) for granted, one can express the first-passage probabihty at 
the level x at time t as 

d °° 

- -T-R{t, x) = Y^ R{2n, x)FToh {Nt = 2n + 1} . (4.10) 

^^ n=0 

Using the fact that the Laplace tranforni {s ^^ i = \nt) of Prob {Nt = n} is equal to 
/"(s) [1 — /(s))]/s, and that the asymptotic form of the generalized persistence probability 
after a fixed number of spin-flips is R{2n,x) ~ [p(x)]", ( 4.1CI|) boils down to a geometric 



summation in Laplace space. A non-trivial pole emerges at s = —{6{x) + 1), where 6{x) 
obeys 

f{-9{x)) = -^. (4.11) 

'p{x) 



It is interesting to remark that ( [4.11D is exactly the equation derived (also within the 



IIA) in [^ for the exponent 9{p) associated with the survival probability of a spin which 
is reborn with probability p each time it flips, providing that we make the correspon- 



dence p = J p{x) . There is however a difference between this extension of the notion of 
persistence and the one we are studying here. In our case, the survival probability p{x) 
imposed every other timestep is not an external parameter, but is directly inherited from 
the dynamics. Nevertheless, one expects that the difference between these two aspects 



should vanish when the survival probability (be it either p or J p{x)) is close to 1, thus 
constituting another plausible argument backing ( [4.9|) ||39|] . 



Now, as 6{x) -^ when a; ^ — 1, the expansion of the Mellin-Laplace transform 
/(s) ^ /(O) + s/'(0) = 1 — s (/) for a small argument s = —9{x) gives 



e{x) ^ {1}-' [1 - ^p{x)j (4.12) 

For the case fz{z) = 6z^^^, where we have an exact expression for p{x), the coefficients 
recombine neatly (using in turn ( |3.18| ) and (^TTUp) to give 



,-1.. . . ..,'l + x 



e{x)^-B-\e,e + i)^-^-j ^e J dx'fMix'), (4.13) 

an equality valid again to lowest order in x — i> — 1. This seems to be the only case where 
this relationship, already found in |^, holds true (thus correcting a slight imprecision 



made in p6[). Otherwise, using (|3.21|) , one has simply an asymptotic proportionality: 



e{x) ~ (1 + x)^ ~ / dx7A/(x'), X -^ -1. (4.14) 



-1 



This proviso made, we believe that, in the limit x — > — 1, the relationship ( [4.14| ) be- 
tween the spectrum of generalized persistence exponents and the limit law of the mean 
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magnetization should hold for any coarsening system, independently of the underlying 
nature (smooth or non-smooth) of the process. Indeed, sites satisfying ( p.l2| ) must at 
least have their local magnetization M > x, and in the long-time limit there is a fraction 
J^dx'fMix') of such sites. Conversely, the decay rate 6{x) (with respect to the loga- 
rithmic timescale i = Int) of the generalized persistence probability R{t,x) should be 
proportional to J^-^^dx'fM{x'), at least when correlations and constraints can be neglected, 
that is for x ^ —1. The simplest possibility (on dimensional grounds) is then that the 
proportionality factor between 6{x) and j!^idx'fM{x') be given by the usual persistence 
exponent 6. However, as testified by ( [4.12[ ), other constants may show up when 6 does 
not parametrize completely the distribution of intervals between zero-crossings (i.e. when 
the latter is not exactly equal to /(/) = 6e~^\ which implies in particular that (/) and 
1/9 are different), which is the generic situation for smooth processes. Note also that, for 



the Levy-model studied in p^, there is no other scale than the one given by 6 (when the 
Levy laws have an infinite mean, i.e. for an exponent ^ < 1), and — thanks to the Sparre 
Andersen theorem (see, e.g., the appendix of [^]) — the equality 6{x) = 9 J^^dx' fuix') 
holds exactly for any level x. 



4.3 Comparison with spectra measured for the diffusion equa- 
tion 

We now turn to a comparison of our analytical results with numerical measurements of 
9{x) conducted for the diffusion equation (Figure H), in one and two space dimensions. 

The behaviors predicted by (|4.^ )-( ^yT^ ) are confirmed, even though they have been 



derived within the IIA. At the lower edge of the spectrum, a power-law rise 6{x) ~ (1+x)^ 
is compatible with the data (despite the obvious numerical imprecisions in this region). 
At the upper edge of the spectrum, where we have a much smaller statistical uncertainty, 
6{x) approaches its limiting value 6 ^ 0.121(1) with a vanishing slope. A numerical fit 
of the data of the form 6{x) — ^ ~ (1 — x)'^ would give an exponent /i ^ 2.2 ± 0.1, which 
is in fair agreement with our prediction ([4.8|) . We have also represented on Figure ^ the 
spectrum of generalized persistence exponents obtained by direct simulation of a temporal 
process for which the laws of successive spin-flip ratios are taken to be uncorrelated and 
distributed according to a pure powerlaw fz{z) = 9z^~^, with an exponent 9 = 0.121. For 
such a zero-dimensional process, clean powerlaws extending over fifteen decades of time 
can be obtained without much effort, and the numerical accuracy on the 6{x) spectrum 
is very high. Let us note that in this case the probability density function fz{z) does not 
vanish near z = 1 and, browsing through the chain of arguments leading to ([4.8|), the 
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Figure 2: 6{x) as measured numerically for the Id diffusion equation (upper curve, full 
line, system size 2.10^), and for a process with uncorrelated ratios of flip-times, powerlaw 
distributed according to fz{z) = 9z^^^, with 9 = 0.121 (lower curve, dashed line). 500 
values for x are considered in both cases. 



correct prediction near x — *> 1 is that 6{x) reaches its limiting value with a finite slope. 
This is indeed observed in our simulations. One also remarks that, even within this 
"doubly" simplified model, the overall agreement with the real spectra is rather good. 

To quantify better the effect of temporal correlations, we have also determined nu- 
merically the 6{x) spectrum for a process with a distribution of (uncorrelated) intervals 
given by the full law (|2.1| ) . We have first followed |jl3| , Taylor-expanding (|2.1j ) in powers 
of 1/s. Inversion term by term gives the expansion of /(/) near / = 0, and the resulting 
series is resummed by constructing rational Fade approximants of the form IP{1'^)/Q{P), 
P and Q being polynomials of respective degree N and A^ + 1. For the 2d diffusion 
equation, we have extended the order of the Fade approximants to A^ = 18. The spin-flip 
ratios are generated using a rejection method, and the 6{x) can then be obtained. A 
comparison with the real spectra computed for the 2d diffusion equation is presented on 
Figure ^. 

On the whole range of the x-values, the agreement is excellent. Of course, beyond 
a certain range of times, the Fade approximants are not able anymore to mimic the 
exponential decay of the probability distribution function f{l) (see inset of Figure |^), 
which is the reason why the curve corresponding to the Fade method is systematically 
overshooting the real one. 

To circumvent this limitation, we have eventually parametrized the probability dis- 
tribution function /(/) as /(/) = {lU{l)/V{l))exp (—91), U,V being polynomials in / of 
respective order A^ and A^ + 1. The coefficients of these polynomials are determined to 
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Figure 3: 6{x) as measured numerically for the 2d diffusion equation (lower curve, full 
line, system size 11585^), and for a process with uncorrelated ratios of flip-times, obtained 
via a rational Fade approximant of ( p.l| ) (upper curve, dashed line). 500 values for x 
are considered in both cases. The inset shows the local slope of the generated /(/) in a 
log-log scale, the horizontal dashed line representing the theoretical value ^ha ~ 0.1862. 
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Figure 4: 6{x) as measured numerically for the 2d diffusion equation (upper curve, full 
line, system size 11585^), and for a process (lower curve, dashed line) with uncorre- 
lated ratios of flip-times, of which the probability distribution function fz{z) = f{l = 
In (l/z))|d//dz| is tailored to match both the small-/ and large-/ behaviors of the law /(/) 
as encoded in (|2.1|). The inset shows the local slope of the generated /(/) in a log- log 
scale, the horizontal dashed line representing the theoretical value 6'nA ~ 0.1862. 
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match faithfully both the small-/ (/(/) = ail + • • •) and large-/ behavior (/(/) ~ aoe^^') 
of the law of the intervals. The implementation of this procedure (with A^ = 4) for the 
2d diffusion equation leads to a discrepancy (Figure 0) on the whole spectrum of 6{x) 
exponents which does not exceed the difference between — and is of the same sign as 
— the numerically measured value 9{= 0{1)) ~ 0.1875 and the value 6iia = 0.1862 . . . 
predicted by the IIA. 

5 Conclusion 

We have shown in this work how it is possible to gain a better understanding of persistence 
properties for processes having a finite density of zero-crossings. The analytical results we 
have obtained, even though they have been derived within an approximate scheme, are 
in excellent agreement with numerical simulations conducted for the diffusion equation, 
which is considered as one of the simplest — yet non-trivial — model of coarsening. 
Furthermore, unexpected connections with relationships obtained in another class of 
stochastic processes have also emerged, which indicates that a certain universality of 
first-passage properties may exist in those systems, at least on a formal level. Let us also 
emphasize that, if one forgets the relationship ( p.l[ ) between the distribution function of 
intervals and the correlation function, it becomes possible [^] to account for all types of 
6{x) spectra encountered in coarsening systems. 

Finally, outside the realm of persistence and of its siblings, the method we have pre- 
sented in Section 3 should have a wide range of applicability for multiplicative stochastic 
processes. For instance, the calculation we have presented in Section 3 for generalized 
persistence properties after a fixed number of spin-flips can be rephrased in terms of 
the exact solution of a random fragmentation problem, the latter topic being of ongoing 
interest |^, ^, ^. It is also striking to note that first-passage exponents calculated 



recently for reaction-diffusion problems in the presence of disorder with the help of an 
asymptotically exact renormalization group also obey hypergeometric equations similar 

to o n. 



References 

[1] Majumdar S N 1999 Current Science (India) 77 370 {cond-mat/9907407) 
[2] Derrida B, Bray A J and Godreche C 1994 J. Phys. A: Math Gen. 27 L357 



18 



[3] Marcos-Martin M, Beysens D, Bouchaud J-P, Godreche C and Yekutieli I 1995 
Physica A214 396 

[4] Yurke B, Pargellis A N, Majumdar S N and Sire C 1997 Phys. Rev. E 56 R40 

[5] Derrida B 1997 Phys. Rev. E55 3705 

[6] Cueille S and Sire C 1997 J. Phys. A: Math Gen. 30 L791 

[7] Cueille S and Sire C 1998 Eur. Phys. J. 5 7 111 

[8] Drouffe J M and Godreche C 1998 J. Phys. A: Math Gen. 31 9801 

[9] Bray A J, Derrida B and Godreche C 1994 Europhys. Lett. 27 175 

[10] Derrida B, Hakim V and Pasquier V 1995 Phys. Rev. Lett. 75 751; 1996 J. Stat. 
Phys. 85 763 

[11] Majumdar S N and Sire C 1996 Phys. Rev. Lett. 77 1420 

[12] Majumdar S N, Sire C, Bray A J and Cornell S J 1996 Phys. Rev. Lett. 77 2867; 
Derrida B, Hakim V and Zeitak R 1996 Phys. Rev. Lett. 77 2871 

[13] Derrida B, Hakim V and Zeitak R 1996 Phys. Rev. Lett. 77 2871 

[14] Oerding K, Cornell S J, and Bray A J 1997 Phys. Rev. E56 R25 

[15] Sire C, Majumdar S N and Rudinger A 2000 Phys. Rev. E61 1258 

[16] Ohta T, Jasnow D and Kawasaki K 1982 Phys. Rev. Lett. 49 1223 

[17] Mazenko G F 1989 Phys. Rev. Lett. 63 1605 

[18] Bray A J 1994 Adv. Phys. 43 357 

[19] Majumdar S N, Bray A J, Cornell S J and Sire C 1996 Phys. Rev. Lett. 77 3704 

[20] Dornic I and Godreche C 1998 J. Phys. A: Math Gen. 31 5413 

[21] Majumdar S N and Bray A J 1998 Phys. Rev. Lett. 81 2226 

[22] Newman T J and Toroczkai Z 1998 Phys. Rev. E58 R2685 

[23] Krapivsky P L and Ben-Naim E 1997 Phys. Rev. E56 3788; Ben-Naim E and 
Krapivsky P L 1998 J. Stat. Phys. 93 583 

[24] Bauer M, Godreche C and Luck J M 1999 J. Stat. Phys. 96 963 

19 



[25] Slepian D 1962 Bell Syst. Tech. J. 41 463 

[26] Baldassarri A, Chate H, Dornic I and Lemaitre A, "Interface dynamics and gener- 
alized persistence in phase ordering systems", preprint. 

[27] Baldassarri A, Boucliaud J P, Dornic I and Godreche C 1998 Phys. Rev. E59 R20 

[28] Kesten H 1973 Acta Mathematica 131 208 

[29] Derrida B and Hilhorst H J 1983 J. Phys. A: Math Gen. 16 2641 

[30] de Callan C, Petritis D, Luck J M and Nieuwenhuizen Th M 1985 J. Phys. A: Math 
Gen. 18 501 

[31] Krug J, Kallabis H, Majumdar S N, Cornell S J, Bray A J and Sire C 1997 Phys. 
Rev. E56 2702 

[32] Kallabis H and Krug J 1999 Europhys. Lett. 45 20 

[33] Dhar A and Majumdar S N 1999 Phys. Rev. E59 6413 

[34] Barbe A 1992 IEEE Trans. Inform. Theory 38 814 

[35] Actually, the persistence probability is slightly different from the distribution func- 
tion of intervals, since in the former, one does not condition on having zeros at both 
endpoints of the interval. But the large-time behaviour (oc exp {—6i)) of these two 
quantities is the same. 

[36] Ben-Naim E, Frachebourg L and Krapivsky P L 1996 Phys. Rev. E53 3078 

[37] Derrida B and Flyvbjerg H 1987 J. Phys. A: Math Gen. 20 5273 

[38] Frachebourg L, Ispolatov I, and Krapivsky P L 1995 Phys. Rev. E52 2540 

[39] Indeed, while this article was being accepted for publication, two very recent works 
by T. W. Burkhardt and G. De Smedt, C. Godreche, J.-M. Luck (preprint cond- 
mat/0009001) dealing with another smooth Gaussian process for which our method 
is relevant have been brought to our attention. These authors have calculated exactly 
(independently and by completely different methods) the partial survival exponent 
6{p) = \ — ^ arcsin | for a randomly accelerated particle (i.e of which the position Xt 
obeys cPxt/dt"^ = rjt, rjt being a Gaussian white noise). This stochastic process (see 
[lU for a list of references on this problem), being the integral of Brownian motion, 
has a finite density of zero-crossings 1/ (/) = -\/3/(27r). The fact that our equation 



20 



( [4.12 ) does match their result 9{p) ~ (-\/3/27r)(l— p) in the p -^ 1-hmit (identifying 



p with Jp{x) in this weakly-constrained regime, as we have argued before) gives 
us further credence to the general validity of the results we have obtained in this 
section. 

[40] Krapivsky P L and Majumdar S N (preprint cond-mat/00061lT} 

[41] Monthus C and Le Doussal P 1999 Phys. Rev. E60 1212 



21 



